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Abstract — In this paper, we contribute to the theoretical 
understanding, the design, and the performance evaluation of 
multi-source multi-relay network-coded cooperative diversity 
protocols. These protocols are useful to counteract the spectral 
inefficiency of repetition-based cooperation. We provide a gen- 
eral analytical framework for analysis and design of wireless 
networks using the Demodulate-and-Forward (DemF) protocol 
with binary Network Coding (NC) at the relays and Cooperative 
Maximal Ratio Combining (C-MRC) at the destination. Our 
system model encompasses an arbitrary number of relays which 
offer two cooperation levels: i) full-cooperative relays, which 
postpone the transmission of their own data frames to help the 
transmission of the sources via DemF relaying and binary NC; 
and ii) partial-cooperative relays, which exploit NC to transmit 
their own data frames along with the packets received from the 
sources. The relays can apply NC on different subsets of sources, 
which is shown to provide the sources with unequal diversity 
orders. Guidelines to choose the packets to be combined, i.e., 
the network code, to achieve the desired diversity order are 
given. Our study shows that partial-cooperative relays provide 
no contribution to the diversity order of the sources. Theoretical 
findings and design guidelines are validated through extensive 
Monte Carlo simulations. 

Index Terms — Cooperative Diversity, Relaying, Network Cod- 
ing, Cooperative Maximal Ratio Combining, Performance Analy- 
sis, Diversity Analysis, Unequal Error Protection Network Codes. 



I. Introduction 

RELAYING and distributed cooperation have recently 
emerged as potential candidate technologies for many 
future wireless applications and standards [1]. Fundamental 
design objective of these systems is to maximally protect 
the transmission of some network nodes (usually known as 
sources) while minimizing: i) the extra bandwidth demanded 
to accomplish this protection; and ii) the resources of the 
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network nodes (usually known as relays) willing to help 
the sources [2]. Early transmission protocols for cooperative 
diversity were mainly based on the repetition coding principle 
with transmissions over orthogonal channels [3], [4]. The main 
limitation of these protocols is that the diversity gain comes 
at the cost of low spectral efficiency. Motivated by these 
considerations, during the last few years many solutions have 
been proposed to overcome, at least in part, the throughput 
reduction of repetition-based orthogonal relaying protocols. 
Some examples are distributed space-time coded protocols [5], 
non-orthogonal protocols [6], successive relaying protocols 
[7], [8], shifted successive relaying protocols [9], two-way 
relay protocols [10], and cognitive cooperation [11], Advan- 
tages and disadvantage of these cooperative diversity protocols 
are discussed in [2]. 

More recently, a new family of cooperative diversity pro- 
tocols has been introduced to overcome the throughput limi- 
tations of repetition-based protocols, while still keeping their 
affordable implementation complexity. They are today known 
as network-coded cooperative diversity protocols, as they 
exploit the emerging concept of Network Coding (NC) [12] for 
a better diversity and throughput tradeoff [13]. Some examples 
where the achievable diversity of network-coded cooperative 
diversity protocols is studied are [14]— [19] for erasure channel 
models and [20]-[23] for error channel models. These papers 
have all shown that NC can be especially beneficial to improve 
the spectral efficiency of multi-user cooperative networks, 
where the available relays must simultaneously serve many 
independent sources. For example: i) in [14] and [16], it 
is shown that NC-aided cooperative protocols are useful to 
overcome the accurate frequency and timing synchronization 
requirements of distributed space-time coded protocols; ii) in 
[17], it is shown that network-coded cooperation provides a 
better diversity-multiplexing tradeoff than distributed space- 
time coded, repetition-based, and selection relaying protocols; 
iii) in [20] and [22], it is shown that NC can reduce the number 
of channel uses from 2N$ to TVs + 1 in multiple-access 
single relay networks with TVs sources, while still achieving 
second-order diversity for every source; and iv) in [23], it is 
shown that, in a multi-source multi-relay network setup, the 
network code can be adequately chosen to provide each source 
a different diversity order, and, thus, a different robustness to 
multipath fading. A comprehensive state-of-the-art survey of 
advantages and disadvantages of these protocols is available 
in [23]. 

Motivated by these potential advantages of network-coded 
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cooperative diversity protocols against repetition-based coop- 
erative diversity protocols, in the present paper we contribute 
to the theoretical understanding and the design of multi-source 
multi-relay cooperative networks with generic binary NC and 
non-identically distributed fading over all the wireless links. In 
our system model, Ns sources are assumed to broadcast, in 
orthogonal time-slots, their data to the available Nr relays 
and to the single destination. The relays first demodulate 
and then perform NC on the estimated data by applying 
a Demodulate-and-Forward (DemF) relaying protocol [24], 
[25]. More specifically, NC is applied to all received packets 
regardless of correct or incorrect data demodulation (error 
channel model). Each relay can apply a different binary encod- 
ing vector to the received packets, thus being able to apply NC 
only on the packets received from a subset of sources. In the 
present paper, it is shown that this degree of freedom enables 
the sources to achieve unequal diversity orders, which can 
be useful for application to heterogenous networks [26] and 
green communications [27]. From our analytical framework, 
guidelines for the design of network codes to allow each 
source to achieve the desired distributed diversity order are 
derived. 

Two classes of relay nodes are considered in the present 
paper: full-cooperative and partial-cooperative relays. The 
relays in the first class are willing to help the sources by de- 
laying the transmission of their own data frames. These relays 
demodulate the packets received from the sources and then 
apply binary NC on a subset of received packets. The relays 
in the second class are willing to relaying the packets received 
from the sources if and only if the transmission of their data 
frames is not delayed. To this end, these relays demodulate the 
packets received from the sources and then apply binary NC 
on both a subset of these received packets and the first packet 
available in their own buffer. Thus, full-cooperative relays 
transmit only redundant bits, while partial-cooperative relays 
transmit, on the same channel use, redundant bits and new 
information packets. As a consequence, partial-cooperative 
relays entail no reduction of the network rate. The investigation 
of how these two classes of relays contribute to the diversity 
order of each source is motivated by the recent results in 
[28] and [29], where partial-cooperative relays are shown to 
be useful to avoid dedicated resources to forward only the 
packets of the sources, as well as to reduce the transmission 
delay of their own data frames. However, in [29] the impact 
of partial-cooperative relays on the achievable diversity of the 
sources is not explicitly investigated. Furthermore, our study 
is motivated by [12] and [30], where, to achieve the network 
capacity, all network nodes are assumed to encode the data 
available in their buffers with the incoming data transmitted 
from other nodes. However, [12] and [30] focus their attention 
on the solvability of the NC problem, and, on the other hand, 
diversity is not investigated. Finally, we mention that the data 
frames available in the buffers of the partial-cooperative relays 
can be the data frames received, in proceeding time-slots, 
from other (than the Ns of the cooperative network) sources. 
Understanding how the partial-cooperative relays contribute 
to the diversity order of the Ns sources is important to 
understand whether they can help these sources without being 



dedicated network elements. Finally, the destination combines 
all the packets received from sources and relays by using 
a generalized Cooperative Maximal Ratio Combining (C- 
MRC) receiver [22], [24], which accounts for both the DemF 
relaying protocol and the binary encoding vectors used at the 
relays. For analytical tractability, it is assumed that sources 
and relays use binary modulation, i.e., Binary Phase Shift 
Keying (BPSK). Even though this assumption might seem 
restrictive, it is worth mentioning that BPSK modulation is 
currently used in many wireless standards, such as WiMax 
(Worldwide Interoperability for Microwave Access) and LTE 
(Long Term Evolution) for control channels. Furthermore, 
it is commonly used for first analytical investigations of 
very complex communication systems. The generalization of 
mathematical framework and analytical diversity assessment to 
non-binary modulations and non-binary encoding vectors is 
currently under investigation. Some comments are provided 
throughout the manuscript. Some preliminary Monte Carlo 
simulation results are available in [31]. Furthermore, recent 
simulation results including channel coding are available in 
[32]-[35]. In particular, recent results in [35] have shown that 
channel coding does not contribute to the diversity order of 
the sources if the channel is assumed to be quasi-static over 
all the wireless links. 

The main findings of our analysis can be summarized as 
follows: i) the diversity order of each source is given by the 
separation vector [36] of the distributed network code, and 
linear block codes with Unequal Error Protection (UEP) prop- 
erties [37] can be used as network codes even in the presence 
of demodulation errors at the relays; ii) only full-cooperative 
relays contribute to the diversity order of the sources. On the 
other hand, partial-cooperative relays contribute neither to the 
diversity order nor to the coding gain of the sources, irre- 
spective of the binary encoding vectors; iii) properly designed 
(deterministic) binary network codes provide greater diversity 
orders than random binary network codes, which achieve 
only first-order diversity; iv) the C-MRC receiver provides 
error performance very close to the Maximum-Likelihood 
(ML) optimum receiver, but with reduced signal processing 
complexity; v) the C-MRC receiver with BPSK modulation 
is shown to satisfy the so-called uniform error property [38], 
which greatly simplifies the analysis of the error performance 
and the design of binary network codes with UEP properties; 
and vi) compared with repetition-based cooperative relaying, 
network-coded cooperative diversity protocols are shown to 
achieve a larger range of diversity orders for the same number 
of total channel uses. 

Finally, we emphasize that compared with state-of-the- 
art papers, which investigate the design and the analysis of 
network-coded cooperative diversity protocols under an error 
channel model [20]-[23], the novelty and contributions of the 
present paper are as follows: i) in [20]-[22], the analysis is 
restricted to a single relay node that combines the packets of 
all the sources. In the present paper, we consider an arbitrary 
number of relays, each of them using a different encoding 
vector to selectively applying NC on the packets received 
from the sources according to the desired diversity order. In 
particular, the system model in [20], [22] can provide only 
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second-order diversity, while our system model guarantees a 
larger range of diversity orders for each source. Furthermore, 
two classes of relays are considered, and their role on the 
end-to-end diversity is discussed; and ii) in [23], a different 
and sub-optimal demodulator is investigated, as well as only 
full-cooperative relays are considered. In the present paper, 
we show that the C-MRC receiver not only provides better 
performance than the demodulator in [23], but it also needs 
lower signal processing complexity. 

The remainder of the present paper is organized as follows. 
In Section [Til system model, transmission protocol, and C- 
MRC receiver are introduced. In Section [III] the error per- 
formance of network-coded cooperative diversity protocols 
is studied, and an asymptotically-tight analytical framework 
is developed. In Section IIV1 the diversity order is studied 
and guidelines for network code design are given. Further- 
more, the role played by full- and partial-cooperative relays 
is discussed. In Section [V] numerical results are shown to 
substantiate our analytical derivation and findings. Finally, 
Section [VT] concludes this paper. 

II. System Model and Transmission Protocol 



buffer. Thus, full-cooperative relays transmit only redundant 
data packets (parity bits), while partial-cooperative relays 
transmit redundant and new data packets. Among the Nr 
relays, N^ c and N^ c act as full- and partial-cooperative 
relays, respectively, with N^ c +N^ C — Nr. The sets of full- 
and partial-cooperative relays are denoted by Nr C and Mr C , 
respectively. It is worth noticing that in our system model the 
relays do not listen, in the relaying phase, to the transmission 
of other relays. In other words, relay-to-relay transmissions 
are neglected in our protocol, and, thus, in our performance 
and diversity analysis. This makes our communication protocol 
sub-optimal. In fact, the parity bits transmitted by the relays 
may be exploited by other relays to provide better estimates of 
the data transmitted from the sources. This option is, however, 
beyond the scope of this paper and its analysis is postponed 
to future research. 

be the BPSK-modulated signal 
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A multi-source multi-relay network with N s sources (S t (XOR) operations; ii) g 
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for t = 1,2,..., Ns), Nr relays (R q forq— 1,2,..., Nr), the 1 x Ns binary encoding vector at R q [12], with gs t R q G 



and a single destination (D) is considered. Transmissions of 
sources and relays occur in orthogonal time-slots. In time- 
slot T f , the source St broadcasts its data packet to D and to 
the Nr relays. This transmission (broadcast) phase lasts Ns 
time-slots. The signals received at R q and D are ys t R q — 
\fEs~t h S t R q xs t + n s t R q and y St D = \fE~s~ t h St DXs t + n StD , 
respectively, where: i) xs t = 1 — 2bs t is the BPSK-modulated 
signal transmitted by St', ii) bs t € {0, 1} is the bit emitted 
by St, iii) Eg t is the transmitted energy per bit of St', iv) 
hxY is the fading coefficient from node X to node Y, which 
is a circular symmetric complex Gaussian Random Variable 
(RV) with zero mean and variance o\ Y jl per real dimension 
(Rayleigh fading). Independent but non-identically distributed 
(i.n.i.d.) fading is considered; and v) nxy is the complex 
Additive White Gaussian Noise (AWGN) at the input of node 
Y and related to the transmission from node X. The AWGN is 
independent and identically distributed (i.i.d.) with zero mean 
and variance No/ 2 per real dimension. 

Upon reception of ys t R q m time-slot T t , the relay R q 
applies ML-optimum demodulation: 



is the data available in the buffer of R, 
= if R q G A/£ c , while g RqRq = 1 if 



{0,1}; iii) b Rq 
and iv) g RqRq 

R q G Nr C . Full-cooperative relays are allowed to transmit 
their own data packets in the first available time-slot at the 
end of the cooperation (broadcast-plus-relaying) phase, which 
lasts Ns + Nr time-slots. Thus, the signal received at D 

where Er is the 



V h R q n R q DX Ra 
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where b s ^ 9 ' denotes the estimate of bs t at R q . 

After Ns time-slots, R q takes turn transmitting, in time-slot 
T/v s + g , a data packet to D. This transmission (relaying) phase 
lasts Nr time-slots. Two classes of relays are considered: full- 
cooperative and partial-cooperative relays. Full-cooperative 
relays transmit a linear combination of the data packets 
received from a subset of sources and delay the transmission 
of their own data packets. Partial-cooperative relays transmit a 
linear combination of the data packets received from a subset 
of sources and the first available data packet in their own 



(1) 



IS VR q D 

transmitted energy per bit of R q . 

In our system model, Channel State Information (CSI) is 
available at the receiver but not at the transmitter. Thus, 
a uniform energy allocation scheme at the sources and at 
the relays is assumed. In formulas, Eg t — E m for t = 
1,2,..., N s , E Rq = E m if R q G A/£ c , and E Rq = E m /2 
if R q G Af£ c for q = 1, 2, ... , Nr. This energy allocation 
scheme takes into account that full-cooperative relays must 
split their available energy to help the sources and to transmit 
their own data at the end of the cooperation phase. 



A. Diversity Combining at the Destination 

Upon reception of all signals ys t D an d Vr d m time- 
slots T t and Tn s+q , respectively, D jointly estimates the data 
transmitted from the sources and the N^ c relays. The data of 
the Nft C relays is transmitted at the end of the cooperation 
phase, i.e., from time-slot Tjv g +jv R +i. Since the data of the 
full-cooperative relays is independent of the data transmitted 
by the sources and by the partial-cooperative relays, it can be 
estimated individually, without loss of optimality, by using a 
demodulator similar to ([T). Demodulation of full-cooperative 
relays is not considered in the present paper, as it can be 
found in textbooks [39]. Let and be 1 x Ns and 



1 x AT 



(PC) 



,, vectors whose entries are b^ for £ = 1,2,..., Ns 



and 6^f ) for R q G JV R 



(PC) 



and q = 1,2,..., A^, respectively. 
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More specifically, and denote the estimates at D of 
the information bits, h$ and b R , actually transmitted from the 
Ns sources and the N R C relays. The vectors bg and h R have 



(D) 



similar definition as b<y and h y R ' , respectively. By assuming 
an ML-optimum diversity combiner, they can be computed as: 



argmax 

6 St e{o,i}, t=l,2 
6 



{a< ml » (b s ,b fl )} 



e{o,i}, R q &Afg c} 



(2) 

where bs and bp are 1 x Ns and 1 x vectors 
whose entries are bg t and respectively, and A^ ML ) (•, •) 
is shown in <(3j at the top of this page, where f^ N< 

gs 1 R q b Sl ®gs 2 R q b S2 <z 

gSxRobSi © 9S 2 R q bs 2 



■®gs Ns R q bs Ns 



Rq 

& (NC) 



gs* 



gR q R q bR q , u R ^ 
R q bs Nv © gR q R q bR 



is 



the network-coded bits at R„ in the absence of demod- 



ulation errors, and P R NC ^ 



Pr 



lN °V&£ C) } is the 



probability that R q forwards a wrong network-coded bit to 
D, which is shown in © at the top of this page, where: 



i) ls t R q 



\hs t R q \ (Eg t /No) is an exponential RV with 



parameter r ys t R q — cr' StR St I No) and probability den- 
sity function f 7StRg (f) = (l/7s t ij,)exp(-^/7 StR J; ii) 
Q (x) = (l /V2w) J x +0 ° exp (-t 2 /2) dt is the Q-function; iii) 

(a) 

= follows because the data of the partial-cooperative relays 

(b) 

are not affected by demodulation errors; iv) = follows from 

0) 

[23, Proposition 1]; and v) « holds for high Signal-to-Noise- 
Ratio (SNR), where ]Jr=t+i I 1 " ^9s r R g Q (V^fsZ^)] -> 
1. 

Even though ML-optimum, the diversity combiner in (0 
is computationally complex. Thus, similar to [22] and [24], 
low-complexity diversity combiners are needed. We consider 
a low-complexity diversity combiner that is best known as 
C-MRC [24], and provide the main steps of its derivation for 
the specific transmission protocol under analysis. According to 
[24, Section II-B], the C-MRC can be derived by regarding 



P R NC) in © as the error probability of an equivalent point- 

u(NC) 



to-point, rather than multipoint-to-point, channel with b R 
at its input and b^ C ^ at its output. In formulas, P R ~ 

E£i gs t R q Q {yJ^ls t R q ) ~ Q (v /2 7e 9 flj, where leqRq is 




10 15 20 25 30 



Fig. 1. Mean and standard deviation of the error function in (|6). Setup: 
i) 9S t R = 1 for t = 1, 2, . . . , Ng; and ii) channel fading is i.i.d. with 
°S t R, = lfor* = l,2,...,AT s . 



the SNR of the point-to-point equivalent channel at R q . 'jeqR 
can be obtained as: 



N s 

E 

t=i 



gs t R q Q (y/2lS t R q ) 



(a) 

> max 

t=l,2,...,Ns 

(b) 



Q 



{Q(^9s! Rq ^ t R q )} 
{9StR q lS t R q } 



(5) 



2 min 

t=l,2,...,Ar s 



(a) 



where: i) > holds because the Q-function is a non-negative 

(b) 

function; and ii) = holds because the Q-function is monoton- 
ically decreasing. Thus, from P R NC ^ ~ Q (^2jeqR q ), we ob- 
tain -f eqRq = min t= i i2i ...,Ar s {gs~tR q Js t R q }, which is an expo- 
nential RV with parameter % qRq = (j2?J x 

gs t R q Js t R q ) 

We note that g^ R = oo if gs t R q = 0. This implies that 
bs t is not network-coded at relay R q , and, thus, the related 
demodulation error probability has no impact in (0. 
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The accuracy of the bound in © is studied in Fig. Q] as a 
function of Ns and SNR. In particular, Fig. Q]shows mean and 
standard deviation of the error function defined as follows: 



Error : 



Ns 

E 

t=i 



9s t R q Q 



S t R„ 



(6) 



Q 



2 min 

4=1,2,. ..,Ns 



> 



Figure Q] clearly shows that §5^ is increasing tight with 
the SNR, and, thus, it is useful for high-SNR analysis, e.g., 
to study the diversity order. Furthermore, we note that, as 
expected, the accuracy of the bound decreases with the number 
of network-coded sources. The error function for Ns = 1 is 
not shown as it would be equal to zero. 

By using the same arguments as in [24, Section II-C], the 
C-MRC with binary NC at the relays is: 



III. Asymptotic End-to-End Error Performance 
of C-MRC 

In this section, the end-to-end error performance of the C- 
MRC in ([8]) is studied. The analytical development comprises 
two steps: i) first, the Average Pairwise Error Probability 
(APEP) is introduced and computed in Section IIII-Al and ii) 
then, the Average Bit Error Probability (ABEP) of the data 
transmitted from the sources and the partial-cooperative relays 
is estimated in Section ITlI-BI and in Section [TlI-Cl respectively. 



A. Computation of the Average Pairwise Error Probability 
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be the hypothesis in (8). Codewords c and c are 1 x (Ns 
Nn) binary vectors, and their entries are c„ and c„ for 
V = 1,2,.. 



argmm 

bs t e{Q,l}, t=l,2,...,N s 

b Rq e{o,i}, R q aM^ c) 



(7) 

with A (C_MRC ) (•, •) shown in © at the top of this page, 
where X Rq = min {j eq R q , lR q o} /jr„d, and j RqD = 
\h R d\ (E Rq /No) is an exponential RV with parameter 
lR q D — a \ q D (E Rq /No). In what follows, we use the 
notation xs t = 1 



2b St and x Rq = 1 - 2b 



It is worth noticing that the C-MRC in (© is different 
from and less computationally complex than the two-step 
demodulator in [23, Eq. (6)]. In Section IVl it is shown that the 
C-MRC in (O also outperforms the two-step demodulator in 
[23, Eq. (6)]. Furthermore, the C-MRC in © is similar to the 
diversity combiner in [22, Eq. (4)]. The difference is that ([8} 
accounts for arbitrary binary encoding vectors, and, thus, each 
relay can apply NC on a different subset of packets received 
from the sources. 



p aim op 

,N S + Nr. The APEP is the probability of 
j^(C-MRC) hiij | detecting c in lieu of c under the assumption that they are the 
only two codewords possibly being transmitted. In formulas, 
APEP(c^c) = E h {PEP(c^c|h)}, where: i) E x {■} 
is the expectation operator computed over the RV X; ii) 
h is a short-hand that denotes all channel gains; and iii) 
PEP (c — > c| h) is the PEP conditioned upon channel fading. 
From (|8]l, PEP (c — > c| h) is given in © at the top of this 
page, where: i) A< C - MRC ) (c,c) = A< C - MRC ) (b s ,b fl ) - 
A (c-mrc) (bsj bfl) . a) g m> with m = 2 o Mo + 2 i m + _ _ + 

2 N «- 1 f i Na ^ 1 and G {0,1} for q = 1, 2, . . . , N R , 

denotes the joint event that the relays with index q such that 
fj,q-i = transmit correct bits, i.e., 6^ 1C ' ) = b^ C \ and 
the relays with index q such that fj, q -\ = 1 transmit wrong 



(NC) 
R q 



bits, i.e., S^ C ' 



7^ b^ R C ^ ■ In formulas, £ m = £m k ^ U £\ 

(a) 



nok) 



shown in dTOb at the top of this page; iii) — follows from 
the total probability theorem by taking into account the 2 Na 
independent demodulation events £ m at the N R relays; iv) 
j = \J — 1 is the imaginary unit; v) n is a short-hand to denote 
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Pv{S m \h}={ 



Nr , ,..„, x N R / N s 

n (i - < c) ) « n i - e gs tRq Q (V2t^) 

g=l v y g=l V {=1 

n< « n E55 t « 3 Q(v^^) 

g=l 9 =1 \t=l / 

Nr / N s \ N R /N s 

~ II U- T, 9s t R,Q (y/2-ystRj X II T,9S t R q Q(y/2'ys t R q , 

9=1 \ t=l / 9=1 \t=l 

Mg-1=0 /Jg-l=l 



if m = 

if to = - 1 

if < to < 2 Nr - 1 



(11) 



a (c-mrc) (c> g | h ) = g (4 7StD 4 t + 4V^d St Re {r& D }) + £ A*, (^zjd^ + 4^5d Rq Re [ V * RqD } 



9=1 



(12) 



A^ a( c-mrc) (s|c,c,£ m ,h) = JJ [cxp (-47 St _Dd St s) exp (47s t Drfs t s )] 

t=i 

JJ [exp (-4A i?5 7 flijjD d^ ) s) exp ^A^jj^rf 



(13) 



2 s 2 



9=1 



(6) 

all AWGNs at D; and vi) = follows from the inverse Laplace 
transform [40, Eq. (5)] with 5 defined in [40, Sec. 2]. 

From © and the independence of the demodulation out- 
comes at the relays, Pr { £ m | h} in ( TTOb is explicitly given 
in ( fTTT i shown at the top of this page. From <[8J and 
some algebra, A' C_MRC ' (c, c| £ m , h) in (O can be sim- 
plified as shown in (fl2l) at the top of this page, where: i) 
Re{-} is the real part operator; ii) Z (•) denotes the phase 
of a complex number; iii) (•)* is the complex conjugate 
operator; iv) r]* tD = (n* StD Ah StD ) / VNo and rf RqD = 
^h RqD ]/^W Q - v) d St = b St - b St G {-1,0,1} 



[ n R q D 
and d Rq = 

= 



l(NC) 
l(NC) 



«4 NC) 

JNC) 



;(NC) ,(NC) 



'Ra 



'R„ 



G {-1,0,1}; and vi) ef£^ 
G {0,1} if =0and 
G {-1,0} if ^9-1 = 1. 



Rq 

From (fT2l . A^ a <c-mrc) ( s| c, c, £ TO , h) = 
E n {exp(-sA( c - MRC ) (c,c|£ m ,h))} in © reduces 
to ( fl3l l, shown at the top of this page, by using the identity 
E n {exp (-KRe{f]* XY })} = exp (K 2 /4) [40, Eq. (19)]. 

Finally, by inserting ( PT3l in (O, the APEP can be re-written 
as shown in ( TBI at the top of the next page, where 7sj? 
and "f RD are short-hands to denote all channels from the 
Ns sources to the N R relays, and from the N R relays to 
D, respectively. For high-SNR, Ft (■, ■) and Q (s| c, c, £ m ) — 
Y[^=i Gq ( s l c ) c, £ m ) can be computed in closed-form as 
shown in Appendix J] The final result is: 

Ji(s|c,c) w [47 StD s (1 - s] 



- bs.-bs 



(15) 



and G q {-\-,-, •) and I (•; •) are given in (TToT l and (fTTI i. respec- 

(o) 

tively, at the top of the next page, and = follows form [39, 
Eq. (5A.9)]. 

From (fT~5b — (fT7T>. A^a( c - mrc ) ( 'I ') i n 02} ca n be written as 



shown in (fT~8T > at the top of the next page. From (fT3Tl-(fT7Ti, 
£ m (c,c) = fajy 1 f g ^^ M A (c-mrc) (s\c,c,£ m ) s~ 1 ds 
in ( fl4l i can be computed by using either the residues theorem 
via the computation of the positive poles of .M A (c-mrc) ( -| •) 
[40, Eq. (6)] or the Gauss-Chebyshev quadrature rule with 5 
being chosen equal to one-half of the smallest real part of the 
non-negative poles of .Ma(c-mrc) (-| •) [39, Sec. 9B.2], [40, 
Eq. (10)]. The residues theorem is to be preferred for simple 
network topologies (with one or two relays) and specific 
network codes that result in expressions of A^a< c - mrc ) ('I ') 
with simple poles. In this case, the integral expression of 
!(•;•) in (TTTb is more convenient to be used compared to 
its closed-form expression. The reason is that closed-form 
expressions of C m (•, •) can be obtained by first applying the 
residues theorem and then computing the integral in (flTt . 
An example is given in Appendix [HI where it is shown that 
our framework coincides with [22, Eq. (32)] for single-relay 
network-coded cooperative networks with encoding vector 
g Rl = [1,1,..., 1]. Another example is available in [41, 
Table I] for one-source two-relay networks. On the other 
hand, the Gauss-Chebyshev quadrature rule is an efficient 
single-integral numerical method that can be used for generic 
network topologies and binary encoding vectors. Its accuracy 
and convergence speed depend on 5, which, from [39, Sec. 
9B.2], [40, Eq. (10)] and by direct inspection of (Q~8]l, can be 
chosen as follows: 



1/2 
1/2 



if 
if 



TO = 



to 7^ and 



5= < 



if to ^ and 



Nr 

E 

9=1 

A»<j-l = 1 

Nr 

E 

9=1 



■ (NC) _ ^(NC) 



'Ra 



Ra 



;(NC) ,(NC) 



Rq 



Rq 



= 

^0 

(19) 
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APEP (C -> C) = 2 £ 1 ^A(C-MHC)(^|c.C, £m ) dg 

m— ^ 



X A (c-mro (s| c, c,£ m ) = E h {A^ A (c-mrc) (s| c, c, £ m , h) Pr{£ m | h}} = J] J" t (s| c, c) x Q (s\ c,c,£„ 
T t (s| c,c) = E 7StD {exp (-47s t _o4 t s ) ex P ( 4 7s t -o4 t s2 ) } 
g (a| c, c, £ m ) = E {7sjj!7flD} | n [ ex P (-4\ RqlRqD d { ™ } s^j exp ^Al^^d^s' 



(14) 



Pr{£ m |7sfl} 



G q (s|c, c,£ m ) 



(47e?fl,s) 1 + (47r 3 ds(1 - s)) 1 -(47egflJ T(s;/x g _i = 0) 

,_1 k(NC) t(NC)| 

(4 7e?Rg ) ^(s^^K ~S I 



j(NO) .(NO) 



if fX q -l = 

if /ig_l = 1 

(16) 



4 z-^/ 2 
X(s;/^_i) = - / l ; _> 



1 



Wo \sin (6») 



+ 4(-l) 



1 (a) I ( 2S ^ 
,s) ,M 



l-V(l + 4s) 1 



-(2,s)- 



if = and s > - (1/4) 



1- V(l-4s)" 



if Hq-\ = 1 and s < 1/4 



(17) 



JVs 



X a( c-mrc) ( s| c, c, £ m ) ps [47 StD s (1 - s)] 



-IH-KI x 



r(NC) ,(NC) 



t=l 



9=1, AI g _l = l 



X II [( 4 7egfl,s) +(4 7i? , ?DS (l- S )) -(47 e?ii J X(s;m 9 -i=0) 



r(NC) ,(NC) 



(18) 



B. Computation of the Average Bit Error Probability of the 
Sources 

From (fl4t and £[8), the ABEP of the sources can be com- 
puted by applying the Union-Bound method [39, Sec. 13.1.3]. 
Let b = [b S) b fl ] and b = [b s , b H ] be 1 x (Ns + N ( R PC) j 
binary vectors defined according to Section IH-AI Then, by 
assuming equiprobable transmitted bits, the ABEP of St is: 



ABEP^ B) < 2 



(a) 



££|d St |APEP(c->c) 



= ^6 St APEP(0 ->■ c) 



(20) 



where |ds t | 



6ft - bf- 



£ {0, 1}, which is introduced in 
(fT2l . takes into account that a wrong decoded codeword, i.e., 
c/c, does not necessarily result in a bit error, i.e., bs t 7^ bs t , 
for S t . 

(a) 

The equality in = , which greatly reduces the computational 
complexity, deserves some comments and a proof. First of 

(a) 

all, = holds for communication systems that satisfy the so- 
called uniform error property [38], i.e., the error probability is 
independent of the transmitted codeword, and, thus, without 



loss of generality, the zero codeword, i.e., c = 0, can be 
assumed to be transmitted. As described in [38], the uniform 
error property does not depend only on the linearity of the 
network code, but it also depends on the code, the modulation, 
the demodulator, and the fading channel. The validity of 
the uniform error property for our system model can be 
proved by direct inspection of (fl8t . More specifically, for a 
given pair (c,c), ( fT8l uniquely depends on: i) the number of 
distinct binary digits between c and c; and ii) the position 
where these binary digits are different. The latter property 
originates from the considered arbitrary network topology 
and the i.n.i.d assumption of channel fading. Let Jp (c, c) = 
{p\ \c p - c p \ = c p ffi Cp = 0, p = 1, 2, . . . , Ns + N R } and 
J P (c, c) = {p\ \dp - Cp\ = Cp © Cp = 1, V = 1. 2, ■ • • ,N S + N R } 
be the sets of Afp and Ns + N R — JVp positions where 
c and c are equal and different, respectively. By direct 
inspection and under the assumption of linear network 
codes, it can be verified that J (c, c) = 3(0, cffic) and 
J (c,c) = (0, c © c) for every pair (c, c). Furthermore, 
since the network code is assumed to be linear, c © c still 
belongs to the codebook. This allows us to affirm that the 
communication system under analysis satisfies the uniform 

(a) . 

error property and, thus, = in (120b can be used to compute 
the ABEP of the sources. 
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>R q = CNs+g © (gs 1 R q b Sl © gs 2 R q bs 2 © ■ ■ • © SSiVgR^Sjvg) =1, #g € A/^ PC) , (c) = 1 j (22) 



M a (c-mrc) (s|0,c, £ m ) = X a( c-mrc) (s|c, £ m ) = (E m /No) 



4f(E)+«4 y (c) + card(£^ ok) ) 



*(s|c,£ m ) 



t&r 1 i-vt 1 - 48 )" 



(5) , 

n 



9=1 



(24) 



ii [Hw) 1 +( 4 xfl,4, D ^(i-^)) 1 '(2s)- 1 ^i- x /(i+4,)- 1 ) 



9=1 
M,-1=0 



r(NC) 



C. Computation of the Average Bit Error Probability of the 
Partial-Cooperative Relays 

With arguments similar to (l20l > and using the uniform error 
property, the ABEP of R q € N R C) is: 



ABEPj£ B) < b Rq APEP (0 -> c) 



E 



APEP (0 -> c) 



(21) 



ceo' 



where ilq PC ^ is defined in ( 122b shown at the top of this page, 



and 



N S +N R - 



is the Hamming weight of c. 



More specifically, tt q ' is the set of possible codewords with 
unit Hamming weight and b R = 1 for R q g N R PG \ The 
reason why only codewords with unit Hamming weight are 
considered is that the diversity order of partial-cooperative 
relays is equal to one. Even though this is intuitive from the 
working operation of our transmission protocol, a formal proof 
of the diversity order of partial-cooperative relays is postponed 
to Section IW^Cl 

(PC) 

By construction, the set fl q contains, at the most, the Ns 
codewords = [0, • • • , 0, 1, 0, . . . , 0] for t = 1,2,..., N s 
where only the t-th entry is non-zero, and the codeword 
gC?) = [0, • • • , 0, 1, 0, . . . , 0] where only the (N s + q)-th entry 
is non-zero. It is worth noticing that is always present in 
fl q PC \ while the presence of for t = 1,2,..., Ns depends 
on the codebook, i.e., the network code. To account for this, 
we introduce the indicator variable £W for t — 1,2,..., Ns, 
such that £W = 1 if £(*) is in the codebook and £w = if cy' 
is not in the codebook. Accordingly, (fJTJi can be re-written as 
follows: 



ABEP ( B UB) < 

tin 



N s 

E 



£ (t) APEP ft) 



APEP 



, , N S 

*=i 



(23) 



(a) 

where ps is proved in Appendix [III] by using the residues 
theorem to compute ([14"1 >. 

From (l23l , the ABEP of the partial-cooperative relays can 
be bounded as (47^)^ + (4 7egi jJ _1 < ABEP^ B) < 

ESl ( 4 7s t D) _1 + (47fl a £>) ~ + (47egijJ ~ ■ Thus, two con- 
clusions can be drawn: i) ABEP^ B) > {ky RqD ) 1 for every 
network topology, network codes, and fading channels. Since 
ABEP ~ (4tr 9 _d) is the high-SNR approximation of the 
error probability of a single-hop transmission from R q to D 
[39], it follows that the partial-cooperative relays undergo a 
performance loss, compared to single-hop transmission, which 
is the price they pay for helping the sources; and ii) this 
performance loss decreases for good source-to-relay channels, 
since in this case j~ qR — > 0, while it increases with the 
number of sources that are network-coded at the relay. In 
Section ITV- Al and in Section ITV-D1 we investigate whether this 
performance loss is offset by a better diversity order and/or 
coding gain for the sources. 

IV. Diversity Analysis and Network Code Design 

The objective of this section is threefold: i) to study the 
diversity order of sources and partial-cooperative relays; ii) 
to provide design criteria for binary network codes such that, 
in addition to having a solvable NC problem, the diversity 
order of each source can be assigned individually (unequal 
error robustness to multipath fading); and iii) to summarize 
relevant performance trends that can be inferred from our 
analytical framework. A byproduct of our diversity analysis 
is a simplified expression of the ABEP of the sources. 

A. Achievable Diversity Order 

We introduce the notation Eg t = xs t Em and E R — 
XR q E m . With the equal energy allocation policy described in 
Section IH we have XSt = h XR q = 1 if Rq G N Rq C \ 

and XR q — 1/2 if R q 6 Af R FC \ By assuming c = 
according to the uniform error property in Section IIII-B1 
A^^(c-mrc) ( • I • ) in (fT~8T > simplified to d24l shown at the top of 

this page, where: i) a 2 eqRg = (E?Ji 9s t R q Xst a s?R q ) I fi ) 

wffi (c) = Et=fi bs t is m e Hamming weight of the systematic 
bits of c, i.e., the bits transmitted from the sources; iii) 
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G d (c,t m ) = w y H ' (c) (c)+ w„ '(cj + w^ (c) = (c) + (c) + (c) 



Gd (c) = min £m {G rf (c, E m )} = min fm ju^f (c) + wffl (c) + w£ ,£m ^ (c) j ( = u^f ) (c) + w^f - 1 (c) 



(25) 



ABEP£ B) « ABEP^ V) = (^ m /iV )- sv ' £ £ (27rj)- 1 / * ( s | E ,f m ) s" 1 ^ (27) 

6, fc St = l, «-H(c)=SV< s t) m=0 V *~ : ' 00 / 



(c) 



fr^ - 1 and w^' £m ^ (c) 



EglJi u 1= i ^iT C ' ) ^ me Hamming weights of the parity 
bits of c, i.e., the bits transmitted from the relays forwarding a 
correct and an incorrect network-coded bit, respectively. Also, 



(5) 



(2) = E»=i 



(NC) 



IS 



,(nok) 



the Hamming weights of all parity bits; and iv) card ^£ 

is the cardinality of £m° k \ i-e., the number of relays that 
transmit an incorrect network-coded bit. 

From (l24l i. it follows that the diversity order of the 
APEP linked to .M^cc-mrc) (-| c, £ m ) depends on c and 



and it is equal to Gd{c,£ m ) 



w 



w 



(S) 
H 



H 



(c) +card (£m° k ^ ) [42]. Since the diversity order, 



G { d Z) , of node Z = (s t , R q £ Af { R P q C) } is determined by 
the APEP with the smallest decaying exponent, we have 



G 



[ P = minc,£ m ^G d (c,£ m )\b z = lj. The conditioning 

upon bz = 1 originates from (120} and (EH . which show that 
the APEP contributes to the ABEP of Z if and only if bz = E 



Let u> 



0,4 : 



complement of 



} (c) - card (d nok) ) - 4f ^ (c) 



be the 



(c). Then, G& (•, •) simplifies to 

(a) 



shown at the top of this page, where: i) = follows 



from card 



(b) 



(£) 



H 



(c); ii) 



follows from mi, (c) = to. 



(2) 



(£); 



iii) in =, the minimum is computed only with respect to £ m 

(d) 

for a fixed codeword c; and iv) = follows from the fact that 

_(R,E (n ° k) ) 

m (c) > by definition. Thus, from ( |25T i: 

G^ Z) =mm{G d (c)| fe z = l} 



bz = 1 



(a) 

= mm 



run | to^ (c) 

un{t»« (c)|6z - l} = SV (z) 



4? (c) 

(6) 



} 



(26) 



(a) 



where: i) = originates from the definition of Hamming weight 



of c, i.e., w<n (c) 



w 



(S) 
H 



(c) + wffl (c); and ii) — follows 



(b) 



from [43, Def. 1] with SV' 2 -* denoting the so-called separation 
vector (SV) of Z. Let a codebook C = {£}, SV (Z) is the 
minimum Hamming weight among the codewords of the 
codebook with bz = E 



In summary, from (|26> two main conclusions can be drawn: 
i) the C-MRC in ([8]) provides node Z with a diversity order 
equal to SV' 2 -*; and ii) by properly choosing the codebook 
and, thus, the network code, each node can achieve a different 
diversity order, which provides it with a different robustness to 
multipath fading and flexibility for improved energy efficiency 
[27]. In Section IIV-CI we provide guidelines to design the 
network code such that the sources can achieve the desired 
diversity order. 

B. Simplified Average Bit Error Probability of the Sources 

The analysis of the diversity order in Section IIV-AI allows 
us to simplify the ABEP in d20l ). In fact, for networks with 
a large number of sources and partial-cooperative relays the 
number of codewords can be very large, which results in the 
computation of many APEPs. Since, for high-SNR, the ABEP 
(l20l is dominated by the APEPs with decaying exponent 
SV (St:) , then (O can be simplified to (f27]i shown at the top of 
this page, where only the codewords resulting in a bit error for 
the sources, i.e., b$ t = 1, and providing the smallest decaying 
exponent with the SNR, i.e., SV^ St \ are retained. 

C. Network Code Design for Cooperative Diversity 

In Section IIV-AI it is proved that sources and partial- 
cooperative relays achieve a diversity order equal to the 
separation vector of the codebook (network code). This result 
is relevant from the system analysis point of view. In this 
section, we look into two important system design issues: i) to 
construct the network code such that the sources have a desired 
separation vector; and ii) to understand which relays contribute 
to the diversity order of the sources. The first design problem is 
relevant because a bigger separation vector could be assigned 
to sources either subject to deeper fading conditions or having 
less residual energy [26], [27]. The second design problem 
is important because nodes that contribute to the diversity 
order of none of the sources provide a marginal contribution 
to their end-to-end performance, and, thus, they may be 
excluded from the relaying phase to reduce the cooperation 
overhead. In our implementation, this decision is made during 
the initialization phase, before actual data transmission, and, 
thus, these excluded nodes will be allowed to transmit the data 
available in their own buffers in the next available time-slot. 
More specifically, three situations can arise for these nodes: 
i) they will act as sources and will transmit their data to the 
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destination; ii) they will act as sources and will ask other nodes 
to act as relays on their behalf for performance improvement; 
and iii) they will act as relays for other sources for which 
they can effectively contribute to their performance. As far as 
this latter case is concerned, a possible scenario is as follows. 
To better understand, we anticipate that partial-cooperative 
relays do not contribute to the diversity order of the sources, 
while full-cooperative relays do contribute. Accordingly, some 
partial-cooperative relays may agree to act as full-cooperative 
relays for other sources, from which they may have received 
special incentives for cooperation. In other words, the relays 
may decide their level of cooperation depending on the sources 
that ask for cooperation. 

Since the diversity order coincides with the separation vec- 
tor of the codebook, it follows that it is uniquely determined 
by the generator matrix, G^ NC \ of the network code. Note 
that this is not true, in general, for the coding gain [42]. 
In fact, it also depends on the demodulation errors at the 
relays and the fading channels. In particular, G^ NC ^ is an 
(n s + N R PC) ^ x (N s + N R ) matrix such that c = bG (NC) . 
By direct inspection of the transmission protocol in Section 
UT1 G( NC ) can be constructed as the block matrix as follows: 



q(NC) = 



(a) 



f-Ns x Ns 





N 



(PC) 



xN s 



I 



N 



(PC) 



xN 



(PC) 







N 



(PC) 



xN 



(PC) 



r (iV s +< C) )x(jV s +< c > 



N 



(PC) 



xN 



(PC) 



(28) 



where: i) I rlX n is an n x n identity matrix; ii) 0; x „ is an 



) is the transpose operator; iv) 
matrix whose q-th column for 



I x n all-zero matrix; iii) 
G( pc ) is the N s x N R PC 

R q € N R PC ^ is the encoding vector g|J ; and v) G^^ is the 
Ns x N R matrix whose q-th column for R q £ Al R is 
the encoding vector g^J . In (|28| ). it is implicitly assumed that 
partial-cooperative relays transmit in time-slots Tpf s+q for 



1,2,. 



,N 



(PC) 



, while full-cooperative relays transmit in 



time-slots T Ns+q for q = N R PC) + l, N {PC) +2, ...,N R . This 
assumption is retained only for simplicity of writing. In fact, 
the performance does not change with the transmission order 
of the relays if the channels are quasi-static during the whole 
cooperation phase, i.e., for Ns + N R time-slots. By direct 
inspection of (1281 l. we notice that G' NC ^ is in row-echelon 

(a) 

form. The identity in = is obtained by applying elementary 
row operations to G ( - NC ^ in order to get a reduced row-echelon 
form. 

From [36, Lemma 4] and [43, Eq. (4)], it is known that 
the separation vector of G^ NC ' can be directly inferred from 
its related N R FC) x (N s + N R ) parity-check matrix H( NC ). 
From (|28]l, H< NC ) is: 



H (NC) 



'xJV^ 



(29) 

are the encoding vectors of the 



where the rows of (G( FC >) 
full-cooperative relays. 

By direct inspection of H^ NC \ and from [36, Lemma 4] 



and [43, Eq. (4)], the following considerations for the design 
of the network code can be drawn: 

• The diversity order of S t is G ( d St) = SV (St) = nc t + 1, 
where ncj is the least number of columns of H^ NC ' whose 
linear combination yields the t-th column of H^ NC \ In other 
words, St has diversity order SV^'-* if and only if the t- 
th column of H^ NC ^ is linearly dependent on no fewer than 
nc 4 = SV (St) - 1 columns of H( NC > [36, Lemma 4]. This 
provides a clear criterion to choose the encoding vectors of 
the full-cooperative relays, and, thus, the matrix G^ PC ' in 
d28j and d29>. 

• Similar to the sources, the diversity order of R q g A/^ PC ' 
is G^ qS> = SV^'-* = ncjv s + 9 + 1, where ncjy s + g is the 
least number of columns of H' NC ' whose linear combination 
yields the (N s + q)-th column of H (NC) . Since the (N s + 
q)-th column of H( NC ' is zero for q = 1,2,..., N R PC \ it 
follows that nc Ns+q = and, thus, G {Pq) = SV^ 3 = 1 
for all partial-cooperative relays. This confirms the analysis 
in Section ITTFCl 

• From d29| i, it follows that the partial-cooperative relays 
do not contribute to the diversity order of the sources. In fact, 
jj(nc) j s independent of the encoding vectors, i.e., the matrix 
q( pc ) in d28l . Only the full-cooperative relays contribute to 
the diversity order of the sources, and, by properly choosing 
the encoding vectors in G^ FC \ the range of achievable 
diversity orders of S t is 1 < SV (St) < N R FC) + 1 [36, Lemma 
5]. 

• In addition to the analytical derivation, it is interesting 
to provide intuition about the result that partial-cooperative 
relays do not contribute to the diversity order of the sources. 
For simplicity, let us assume that there are no decoding errors 
at the relays (i.e., the SNR of the source-to-relay links goes to 
infinity). As far as diversity is concerned, this is the best-case 
scenario. If the diversity order is equal to one in this case, 
then it will be equal to one for every SNR of the source- 
to-relays links. However, this no longer holds if the diversity 
order is greater than one, since decoding errors at the relays 
may reduce the diversity order. To start with, let us consider the 
simplest case study with a single source and a single relay. The 
relay acts as partial-cooperative, and, thus, it network-codes 
the bit received from the source with its own bit. Accordingly, 
the codewords of the codebook are: "00", "01", "10", "11" 
where the first bit is sent by the source and the second bit is 
the XOR of the bits of source and relay. It is apparent that the 
codebook contains all the possible codewords, and, thus, the 
Hamming distance of the information bits can only be equal to 
one. More in general, if all the relays are partial-cooperative, 
i.e., the network rate is equal to one, the intuitive reason why 
partial-cooperative relays do not contribute to the diversity of 
the sources is that the codebook of the equivalent distributed 
code contains all possible codewords of the universe set. As 
a consequence, the Hamming distance of each information 
bit can only be one. If the network rate is less than one, 
i.e., some relays are full-cooperative and some others are 
partial-cooperative, then providing an intuitive explanation is 
more complicated and the equivalent codebook needs to be 
investigated by direct inspection. Our understanding is that 
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when a partial-cooperative relay network-codes its own bit 
with the bits received from the sources the net effect is that 
the transmitted network-coded bit no longer acts as a parity bit 
because of the presence of the information bit of the partial- 
cooperative relay. On the other hand, the transmitted network- 
coded bit is an actual information bit, and, thus, the Hamming 
distance of the equivalent code is reduced compared with 
the case study where all the relays are full-cooperative. Our 
research is currently devoted to understand whether partial- 
cooperative relays can contribute to the diversity order of 
the sources by increasing the Galois field and/or for some 
modulation orders. We are particularly interested in system 
setups with network rates less than one. In fact, if all relays 
are partial-cooperative the equivalent codebooks may coincide 
with the universe set for many choices of Galois field and 
modulation order. If so, the diversity order of the sources will 
never be greater than one. 

From these considerations, the following design criterion for 
the network code is proposed. Let us consider a network with 

(S ) IFO 

Ns sources requesting diversity order G d and N R full- 
cooperative relays. Partial-cooperative relays can be neglected 
since they do not contribute to Gj . Then, G^- FC \ i.e., the 
encoding vectors of the full-cooperative relays, can be chosen 

{ FC} 

as the generator matrix of a systematic (Ns + N R ,Ns) 
linear block code with separation vector SV { *' — G d . 
Known techniques for the construction of these codes, along 
with the conditions for their existence for a given choice 
of Ns and N R , are available in the literature [36], [44]. 
Examples of such systematic codes are listed in [37, Table I] 
for various network topologies and separation vectors. Among 
the possible choices, an important case study is when each 
source needs the largest separation vector, i.e., 
N R FC ^ + 1 for t = 1, 2, . . . , Ns- This case study corresponds 
to the design of the so-called Maximum Distance Separable 
(MDS) codes [44, Ch. 11], which are known to attain the 
Singleton bound [45]. From [44, Ch. 11, Corollary 7], it 
follows that for Ns > 2, binary modulations, and binary 
network codes the Singleton bound can be attained only for 
single-relay networks. If more relays are available, non-binary 
codes should be used to attain the Singleton bound. For a two- 
source two-relay network with binary NC, this result has been 
proved, using simulations, in [26] by considering all possible 
choices of the encoding vectors. On the other hand, for single- 
relay networks the possibility to attain the Singleton bound 
for arbitrary Galois fields has been proved in [22]. In fact, the 
network code used in [22] is an example of the the so-called 
trivial MDS codes described in [44, Ch. 11, Problem 1] with 
code type (n, n — 1,2) and n = Ns + 1. 

D. Insights From the Analytical Framework 

We close this section by providing some conclusions that 
can be drawn from our analytical framework and can help the 
understanding, the analysis, and the design of network-coded 
cooperative networks. 

CI) From the definition of <j\ qR in d24l >. it follows, as 
expected, that the larger the number of network-coded sources 
at each relay, the worse the performance. This result is in 



agreement with [22] and [23], as well as with intuition. 
Mathematically speaking, this originates from the fact that 
increasing the number of network-coded sources results in 
a larger probability, P R NC ^ in (|4]i, of transmitting an incorrect 

network-coded bit. Thus, the matrix (G^ FC ^) T in (|29| | should 
not only be designed to provide the desired diversity order, 
but its rows should be as sparse as possible in order to 
minimize the error accumulation problem at the relays, and, 
thus, to improve the coding gain [42]. In general, the design 
criterion on the diversity order has a more pronounced impact 
on the end-to-end performance. Furthermore, the closed- 
form expression of cr\ qR in d24i i confirms that the error 
accumulation problem at the relays is inversely proportional to 
the quality of the source-to-relay links, i.e., a SR . This result 
is well-known for DemF-based relaying protocols [24], [25], 
and it can be considered as a sanity check for our analytical 
derivation. 

C2) In Section ITV-Cl it is shown that SV (St) > 1 for 
every choice of (G^ c )) in (|29>. In particular, S t enjoys 
no diversity gain if and only if SV' 5 '-' = 1, which implies 
gy(5't) _ nCt ^i = i=^ nC( = o. This is possible if and 

only if the t-th column of (G^ c ') is all-zero, which, in 
turn, implies that St is network-coded by none of the full- 
cooperative relays. This result is expected and agrees with 
intuition. In fact, if St is network-coded by none of the 
full-cooperative relays this implies that it is only present in 
the direct link, and, thus, the diversity order can only be 
equal to one. This sanity check confirms the consistency of 
our analytical derivation. In addition, and, more importantly, 
it provides an interesting conclusion. It is sufficient that at 
least one entry of all columns of (G^ ') is non-zero 
to guarantee that each source enjoys at least second-order 
diversity, i.e., G d St) = SV (St) > 2. In other words, if each 
source is network-coded by at least one full-cooperative relay, 
then its diversity order is at least equal to two regardless of 
the specific choice of the network code, which can also be 
randomly generated. 

C3) Let us assume that the N R FC ^ full-cooperative relays 
have the same encoding vector and that at least two of its 
entries are non-zero, i.e., at least two sources are network- 
coded. If the encoding vectors are all the same and only one 
of its entries is non-zero, the system reduces to a cooperative 
network without NC. This latter case study is considered in 
C6) below. Under these assumptions, the sources that are 
network-coded have diversity order two, while the sources that 
are not network-coded have diversity order one. In fact, the 

rp J 

matrix (G( FC )) in <|29} has a number of all- one columns 
equal to the sources that are network-coded, and the other 
columns are all-zero. The all-zero columns are responsible for 
the first-order diversity of the sources that are not network- 
coded. On the other hand, the all-one columns, since linearly 
dependent, are responsible for the second-order diversity of 
the network-coded sources. 

( FC} 

C4) Let us consider that the N R full-cooperative relays 
network-code the data received from all the sources. This 
results in having an all-ones matrix (G^ c ') in d29| i. Ac- 
cordingly, its columns are linearly dependent and the diversity 
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order of each source is G { d St) = SV (St) = 2. 

C5) In [22], it is shown that one relay is sufficient for all 
network-coded sources to achieve second-order diversity. The 
diversity analysis in Section HV-CI reveals that this encoding is 
optimum from the diversity point of view. In fact, we have 



= 1, 



proved that 1 < SV (St) < 7Vjf C) + 1. If N R FC) 
then 1 < SV (5t) < 2. In [22], it is proved that G { f t} = 
SV (St) = 2, which is the best achievable diversity with 
one full-cooperative relay. From C3), it follows that multiple 
relays can be useless, from the diversity point of view, if their 
encoding vectors are all the same since diversity orders no 
greater than two can be achieved. In the presence of multiple 
relays, the encoding vectors should be carefully designed in 
order to enjoy the benefits of cooperative diversity. The design 
guidelines provided in Section IIV-CI are useful to this end. 

C6) As mentioned in Section Q] network-coded cooperation 
is a generalization of repetition-based relaying. In particular, 
the former reduces to the latter if there is only one non-zero 
entry in each encoding vector gji q for R q € M R FC \ Thus, 
the frameworks developed in the present paper, e.g., ( f2Tb . 
and the diversity analysis in Section IIV-AI can be applied to 
repetition-based relaying as well. In this case, the t-th column 
of (G^* 7 )) in d29l i has a number of non-zero entries equal 
to the full-cooperative relays that forward the data received 
from S t . Let < N { R St] < N R FC) denote this number of 
relays. Also, since each relay can forward the data of one 
and only one source, it follows that the first N$ columns 
of (G( FC )) are independent of each other and that net = 
N R St) . Thus, G ( f t] = SV (St) = nc t + 1 = + 1. This 

result is in agreement with state-of-the-art diversity analysis 
of repetition-based relaying [24]. However, unlike [24], our 
frameworks provide also an accurate estimate of the coding 
gain. 

C7) In Section IIII-C1 it is proved that the partial- 
cooperative relays undergo a performance loss, compared to 
single-hop transmission, which can be interpreted as the price 
they pay for helping the sources. Furthermore, in Section llV-CI 
it is shown that the partial-cooperative relays do not contribute 
to the diversity order of the sources. Thus, a fundamental ques- 
tion is whether the partial-cooperative relays are useful and 
whether they should be exploited during the relaying phase. To 
answer to this question, it is important to investigate the impact 



that these nodes have on the coding gain of the sources. To this 
end, let us consider two network topologies: 1) the first one 
with Ns sources, N R FC ^ full-cooperative relays, and N R PC ^ 
= partial-cooperative relays; and 2) the second one with Ns 
sources, N 1 ^ ^ full-cooperative relays, and ivjj PC ^ partial- 
cooperative relays. In both networks, the full-cooperative 
relays have the same encoding vectors, while in the second 

f PC) 

network the N R partial-cooperative relays have arbitrary 
encoding vectors. Let C\ and C2 denote the codebooks of the 
two network topologies. In general, C\ 7^ C2. From Section 
IIV-CI it is known that, in both networks, the sources have the 
same diversity order, i.e., Svf t} = SV { 2 St) = SV (St) . From 
d27), we know that for high-SNR the ABEP of S t depends 

(S ) 

only on the codewords Ej e Ci with SVj = wu (c*) = 
w u ^ fii) + w jP (c?) for Z = 1, 2. For these codewords, the 
constraint SVp - w(St) ; — - {S) - - {S) 
„(«) 



SV^ t; implies w£> (c*) = w%> (c|) 



and w^ 1 (cj) = w^p (c|), which stems from the fact that 
both codebooks have, by construction, the same systematic 



bits. In addition, we can write (c*) = w^'*'^' (c*) + 
w [P ' PC) (c?), where w [ £' FC) (cf) and w [P ' PC) (£,") are the 
contributions to the Hamming weight related to full- and 
partial-cooperative relays, respectively. By construction: i) 
w (R,FC) _ w (r,fc) j-g*^ since in both networks the 
full-cooperative relays use the same encoding vectors; and 

(R PC) 

ii) w H ' (cf ) = 0, since in the first network there are 
no partial-cooperative relays. In conclusion, the constraint 



.,(«) IT. 



(£!) 



■^f ) (c|) implies w^'^' (c* 2 ) = 0, which, 



(R,PC) 



(NC) = for R q e 



in turn, implies, for these codewords, b q 
■Mr ■ in other words, in the second network, the partial- 
cooperative relays do not contribute to the Hamming weight 
of the codewords ci, S C2. From ^ ( -| c|, •) in d24"l i. we notice 
that the relays with b q NC ^ — do not contribute to the coding 
gain of the ABEP. Thus, the partial-cooperative relays do not 
contribute to the coding gain of the sources. 
In summary, our analytical study clearly shows that the 
partial-cooperative relays contribute to neither the diversity 
order nor to the coding gain of the sources. In addition, 
these relays undergo a performance loss due to applying 
NC on the data received from the sources. In conclusion, 
our analysis reveals that only full-cooperative relays should 
be used during the relaying phase. This result allows us 
to further simplify the computation of \& ( • | • , •) in 
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and ABEPjF^ in dZTT l. Since the partial-cooperative relays 
contribute to neither the diversity order nor to the coding 
gain, they can be completely neglected from the computation 
of #(-|-,-) and ABEP^ V) . Thus, dH and <27) simplify 
as shown in (130b at the top of the previous page, where 



bs 1 ,bs 2 , - ■■, bs A 



r(NC) 

; ' °Ri ■ 



r(NC) 
U R 2 ' 



ac- 



counts only for sources and full-cooperative relays. In other 
words, in the high-SNR the performance of the sources is 
determined by the sub-network containing only the sources 
and the full-cooperative relays. The partial-cooperative relays 
can be neglected regardless of the binary encoding vectors 
they use. 

Finally, some remarks about the usefulness of partial- 
cooperative relays are needed. In the present paper, we have 
focused our attention on the contribution of these relays to 
diversity order and coding gain. Our results have clearly 
shown that these relays have a negligible contribution to 
these performance metrics. However, other papers have shown 
that partial-cooperative relays have many advantages, which 
include energy efficiency, low transmission delay, as well 
as the possibility to avoid dedicated network elements for 
data relaying [29]. Furthermore, partial-cooperative relays are 
always present in the general NC framework [12], [30]. As a 
consequence, our framework has highlighted important perfor- 
mance tradeoffs for different types of relays. Full-cooperative 
relays seem to be more useful for diversity and end-to-end 
performance. On the other hand, partial-cooperative relays 
seem to be more useful for energy efficiency, higher rate, 
and low transmission delay. How these tradeoffs are affected 
by non-binary modulation and non-binary Galois field is 
currently being investigated by the authors. 

C8) So far, the encoding vectors at the relays have been 
assumed to be fixed a priori, i.e., deterministic NC. This 
implies that the encoding vectors must be agreed by relays 
and destination before transmission. A different option is to 
use random NC, which allows the relays to generate at random 
the encoding vectors. This solution is suitable for distributed 
implementations [46], [47]. However, random binary NC can- 
not guarantee that the sources achieve a given diversity order 
since there is no a priori structure on the network code. In fact, 
in random binary NC the worst case scenario for the diversity 
of source St happens when none of the encoding vectors at 
the relays include St, i.e., when gs t Rq — for R q 6 Njf ^ 
and q — 1,2,..., Nr. In C2), it is shown that in this case the 
sources have diversity order one. Since St is network-coded 
in none of the relays, this implies that its error performance 
is equal to the error probability of a single-hop transmission 
multiplied by the probability that gs t R q — for R q £ J^ R FC ^ 
and q = 1,2, ... , N R , i.e., 1 / 2 N R C \ Thus, for high-SNR, 



ABEP 



(SV, random) 



1/2 



N 



(fc; 



[4Xs t o* St D(E m /N )]- 



As a result, random binary NC allows the sources to achieve 
better end-to-end error performance than single-hop trans- 
mission. Also, the performance gain increases with N R FC \ 
However, it does not help increasing the diversity order. 

Finally, we would like to emphasize that the conclusions 
drawn in the present paper depend in many cases on the 



assumptions of binary modulation and binary NC. The im- 
pact of non-binary modulation and non-binary NC on these 
conclusions is currently being investigated by the authors. 

E. Generalization to Non-Binary Modulation and Non- 
Binary Network Code: Thoughts and Conjectures 

Diversity analysis and related network code design dis- 
cussed in the previous sections are applicable only to the 
system setup with binary modulation and binary network 
codes. In the present paper, these assumptions are retained only 
for analytical tractability and to provide sound proofs about 
the achievable diversity for generic network topologies and 
realistic channel models over all the communications links. 
In this section, we provide some thoughts (or conjectures) 
about the impact of the Galois field on the expected diversity 
order. The departing point of our thoughts moves from the 
outcomes obtained in the binary case. We have shown that if 
a proper demodulator at the destination is used, then classical 
systematic linear block codes can be used as network codes. 
More specifically, the design of diversity-achieving network 
codes is equivalent to the design of linear systematic block 
codes over fully- interleaved point-to-point links (see Section 
IIV-CI and [23]). In other words, if the demodulator is well 
designed to account for demodulation errors at the relays, then 
network codes for distributed diversity can be constructed by 
assuming ideal (error-free) source-to-relay links. In what fol- 
lows, we call this demodulator, assuming that it exists for every 
system setup, genie-aided detector. The two-step demodulator 
in [23] and the C-MRC demodulator in Section III- Al are two 
practical examples of detectors for binary modulation and 
binary network codes. Likewise, the C-MRC demodulator in 
[22] for single-relay networks is another practical example 
under the assumption that modulation order and Galois field 
are the same, as well as that a MDS network code is used 
[44, Ch. 11]. Let us now assume that such demodulator exists, 
regardless of its computational complexity and network CSI 
needed at the destination to counteract the error propagation 
problem. Under these assumptions and if modulation order 
and Galois field are the same, we expect that the results for 
the binary case can be generalized to the non-binary case. 
In particular, we conjecture that non-binary linear systematic 
block codes can be used as diversity-achieving network codes. 
This implies that existence conditions and code constructions 
for equal (e.g., MDS code design) and unequal (based on the 
separation vector) diversity apply. For example, from [44, Ch. 
11, Corollary 7] and [45] this would imply that the Singleton 
bound can be achieved for the design of MDS codes if the 
size, q, of the Galois field of the network code (and, thus, the 
modulation order) satisfies the inequalities q > N$ + 1 and 
q > N R FC) + 1 if N s > 2 and N R FC) > 2. This conjecture 
is similar to the results recently obtained in [15], [17]— [19] 
for the erasure channel model. On the other hand, the analysis 
of the system setup where Galois field and modulation order 
are different is more complicated to be conjectured since it 
would need further assumptions on how modulation and NC 
are mixed together. Formal proofs of these conjectures, as 
well as the design of practical demodulators to obtain these 
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Fig. 2. ABEP of a 2-source 2-relay network. Both relays are full- 
cooperative. Performance comparison of three diversity combiners at the des- 
tination: i) the ML-optimum demodulator in {5}; ii) the two-step demodulator 
in [23, Eq. (6)]; and iii) the C-MRC in (8). Setup: i) channel fading is i.i.d. 
with Oq = 1; and ii) the encoding vectors are g^ = [1,0], gfl 2 = [1,1], 
which yield SV< Sl ^ = 3 and SV^ 2 ) = 2. For clarity, only Monte Carlo 
simulations are shown. 



Fig. 3. ABEP of a 2-source 5— relay network. All relays are full-cooperative. 
The C-MRC in {8) is used at the destination. Setup: i) channel fading is i.i.d. 
with Oq = 1; and ii) the encoding vectors are gR 1 = [1,0], gjj 2 = [1,0], 

gi? 3 = [1.1]. Sfl 4 = [1.1]- SR 5 = [0,1], which yield SV< s i) = 5, 
and SV^ 2 ' = 4. Markers show Monte Carlo simulations, solid lines show 
ABEP^ B) , and dotted lines show ABEP^ V) . 



diversity gains for realistic source-to-relays links and for an 
arbitrary number of sources and relays are currently being 
investigated by the authors. 

V. Numerical and Simulation Results 

In this section, analytical framework and findings are 
validated through Monte Carlo simulations. For simplicity, 
numerical examples for i.i.d. fading channels are described. 
However, the framework has been verified for various channel 
conditions. In all analyzed scenarios, we have obtained a good 
accuracy for high-SNR. An example for i.n.i.d fading channels 
is shown in Fig. [6] The encoding vectors are obtained, for 
a given network topology and SV, from [37]. As far as the 
analytical framework is concerned, the ABEP of the sources, 
ABEP^ V) , is computed using (f30b and the Gauss-Chebyshev 
quadrature rule in [39, Sec. 9B.2], [40, Eq. (10)] with 6 given 
in ( fl9] l. For completeness, ABEF^ 8 ' in ( 1201 is shown as 
well. It is obtained from (f30b by neglecting the condition on 
SV (St) . 

In Fig. 12 a network topology with two sources and two 
relays is studied, and the ABEP of three diversity combiners is 
compared. It can be observed that the C-MRC provides near- 
ML performance with reduced signal processing complexity. 
Furthermore, it provides, with lower implementation complex- 
ity, better performance than the two-step demodulator in [23]. 
In fact, the weighting factors of the C-MRC are simpler to be 
computed. Other network topologies with more sources and 
relays have been studied, and the same performance trend as 
in Fig. |2] has been observed in all analyzed case studies. 

In Fig. [3] and Fig. @1 a network topology with two sources 
and five relays is considered, and the analytical frameworks 
are compared with Monte Carlo simulations. For high-SNR, 
we observe a good accuracy. It is worth mentioning that, 
unlike ABEP^ JB , ABEP^ V ' is not an upper-bound since 



some codewords are not considered in the computation. How- 
ever, ABEP^ V ' asymptotically (high-SNR) overlaps with 
ABEP^ B) , as discussed in Section HSU The ABEP of the 
partial-cooperative relays obtained from (l23l is shown as well. 
The figures confirm the correctness of the diversity analysis in 
SectionUV] More specifically, Fig.|4]shows the diversity order 
degradation caused when two relays in Fig. [3] are no longer 
full-cooperative but become partial-cooperative. The ABEP of 
Fig. |4] coincides with that of the sub-network without partial- 
cooperative relays. In Fig. [5] and Fig. [7] the same analysis is 
conducted for a network topology with three sources and three 
relays. Conclusions similar to Fig. |3]and Fig. |4]can be drawn. 
Furthermore, Fig. [6] shows an example for i.n.i.d. fading by 
assuming the same network topology and network code as 
in Fig. [5] More specifically, the source-to-relay channels are 
assumed to be of better quality than the other channels. We 
observe a very good accuracy of our framework. In addition, 
by comparing Fig. [5] and Fig. [6] we observe, as expected, that 
the ABEP in Fig. [6] is better since error propagation is less 
pronounced for this setup. 

In Fig. [8] the ABEP of repetition-based cooperative relaying 
and random binary NC is shown. A network topology with 
three sources and three relays is investigated. Figure |8]con firms 
that repetition-based cooperative relaying is a special case 
of network-coded cooperation, and that our framework can 
be used to provide good estimates of both diversity order 
and coding gain. Also, Fig. [8] confirms the first-order di- 
versity of random binary NC, as well as the accuracy of 
ABEP^ V ' random) computed in Section USE 

In Fig.|9]and Fig.[l0] the ABEP of network-coded (NC) and 
repetition-based (RepCod) cooperative diversity is compared 
by assuming the same network topology (two sources and five 
relays) and the same number of time-slots. Figure [9] and Fig. 
[TOl show that, by properly choosing the encoding vectors, the 
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Fig. 4. ABEP of a 2-source 5-relay network. Ri and R2 are partial- 
cooperative relays, and R3, R4, and R$ are full-cooperative relays. The C— 
MRC in is used at the destination. Setup: i) channel fading is i.i.d. with 



1; and ii) the encoding vectors are gjj x = [1,0], gR 2 



[1.0], 



gi?, = gi? 4 = [i.i]. §i?5 = [0.1]. which y ield sv( 5 i) = 

SV (S2) = 3. As for the sources, markers show Monte Carlo simulations, 



solid lines show ABEP 



(UB) 



and dotted lines show ABEP 



(SV) 
St 



More 



specifically, ABEPg^ B ^ and ABEP^ V ' are computed by considering only 
the sub-network with sources and full-cooperative relays. As for the partial- 
cooperative relays, markers show Monte Carlo simulations, and dotted lines 
are obtained from 423) . 



-* — C-MRC (S1) 
B — C-MRC (S2) 
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Fig. 5. ABEP of a 3-source 3-relay network. All relays are full-cooperative. 
The C-MRC in (SJ is used at the destination. Setup: i) channel fading is 
i.i.d. with Oq = 1; and ii) the encoding vectors are g Rl = [1,0,1], 
gJ? 2 = [1.1.0], gj? 3 = [1,0,0], which yield SV (Sl) = 4, and SV (S2) = 
Sy( s 3) = 2. Markers show Monte Carlo simulations, solid lines show 



ABEP 



(UB) 
St 



and dotted lines show ABEP 



(SV) 
St ■ 



Fig. 6. ABEP of a 3-source 3-relay network. All relays are full-cooperative. 
The C-MRC in (8) is used at the destination. Setup: i) channel fading is 
i.n.i.d. with erg = i t a 2^ = CT 2^ = ff 2 j and CT 2^ = (xooOo-q) 2 

for t = 1, 2, . . . , Ng, q = 1, 2, . . . , Nr, i.e., the source-to-relay links 
are better than the other links; and ii) the encoding vectors are gn 1 = 
[1,0,1], g R2 = [1,1,0], g Ra = [1,0,0], which yield SV< s i) = 4, and 
SV^ 2 ' = SV' S3 ' = 2. Markers show Monte Carlo simulations, solid lines 



show ABEP 



(UB) 
St 



and dotted lines show ABEP 



(SV) 
St 
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Fig. 7. ABEP of a 3-source 3-relay network. R\ is a partial-cooperative 
relay, and R2 and R3 are full-cooperative relays. The C-MRC in {8} is used 
at the destination. Setup: i) channel fading is i.i.d. with <Tq = 1; and ii) 
the encoding vectors are g Rl = [1,0,1], gjj 2 = [1,1,0], g Ra = [1,0,0], 
which yield SV (Sl) = 3, SV (S2> = 2, and SV (S3> = 1. As for the sources, 



markers show Monte Carlo simulations, solid lines show ABEP 



(UB) 



and dotted lines show ABEP 

d(SV) 



(SV) 



More specifically, ABEP^y B) 



and 



are computed by considering only the sub-network with sources 
and full-cooperative relays. As for the partial-cooperative relay, markers show 
Monte Carlo simulations, and dotted lines are obtained from (23). 



transmission protocol described in Section HJis flexible enough 
to accommodate various performance tradeoffs for the sources. 
For example, the setup denoted by NC-3 outperforms the setup 
denoted by RepCod-3. In fact, the ABEP of S\ is almost 
the same and the ABEP of S2 is much better thanks to the 
higher diversity order. This choice of the encoding vectors is 
convenient when S\ needs high robustness to multipath fading. 
On the other hand, the setup denoted by NC-2 appears to be 



useful when both sources require almost the same ABEP. In 
fact, compared with NC-3 and RepCod-3, the performance 
degradation of Si is tolerable, and the performance gain of 
S2 is significant. Furthermore, NC-2 outperforms RepCod-1, 
and it provides a good alternative to RepCod-2. In fact, the 
big performance gain enjoyed by S2 is obtained with a small 
performance loss for Si. Finally, it is apparent that random 
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Fig. 8. ABEP of a 3-source 3-relay network. The relays use repetition- 
based cooperation. The C-MRC in {8) is used at the destination. Setup: i) 
channel fading is i.i.d. with Oq = 1; and ii) the encoding vectors (only one 
source is relayed) are gR 1 = [1,0,0], gR 2 = [1,0,0], gR 3 = [0,1,0], 
which yield SV (Sl) = 3, SV (Sa) = 2, and SV (S3) = 1. Markers show 



Monte Carlo simulations, solid lines show ABEP 



ABEP 



(SV) 



(UB) 
St > 



and dotted lines show 



The figure shows also the ABEP when the relays generate at 
random the binary encoding vectors: random binary NC. In this case, markers 



show Monte Carlo simulations and dotted lines show ABEP 



Section ITV-Dl 



(SV, random) 
St 



- NC-1 - SV=(2,2) 
NC-2 - SV=(5,4) 

- NC-3 - SV=(6,2) 

- Random NC - SV=(1,1) 

- RepCod-1 - SV=(4,3) 
RepCod-2 - SV=(5,2) 
RepCod-3-SV=(6,1) 




Fig. 9. ABEP of a 2-source 5-relay network (source Si). All relays are 
full-cooperative. The C-MRC in {8) is used at the destination. Channel fading 
is i.i.d. with Oq = 1. Various encoding vectors are considered to compare 
network-coded (NC) and repetition-based (RepCod) cooperative diversity: i) 
NC-1 with g Rl = SR2 = g Rs = g Ri = g fl5 = [1,0]; ii) NC-2 with 
gfl! = SR 2 = [1,0], S r 3 = g R4 = [1,1], gR 5 = [0,1]; iii) NC-3 with 
gi?! = [1,1]. SR 2 = &R 3 = SR 4 = gR 5 = [1,0]; iv) RepCod-1 with 
SR! = SR 2 = SR 3 = [1,0]- SR 4 = SRs = I ' 1 ]' v ) RepCod-2 with 
gJ?! = SR 2 = SR 3 = gfl 4 = I 1 ! !' gfls = [°,1], vi > RepCod-3 with 
gfli = SR 2 = SR 3 = SR 4 = SR 5 = [1,0]; and vii) Random NC with 
binary encoding vectors generated at random. The separation vectors of both 
sources are shown in the legend of the figure. For clarity, only Monte Carlo 
simulations are shown. 



binary NC is, in general, not a good solution, especially in 
the high-SNR, because of the single-order diversity order. 

In addition to the numerical examples depicted in Fig.|2UT0l 
the considerations in Section IT V-DI have been verified through 



Fig. 10. ABEP of a 2-source 5-relay network (source $2). All relays are 
full-cooperative. The C-MRC in (8j is used at the destination. Channel fading 
is i.i.d. with Oq = 1. Various encoding vectors are considered to compare 
network-coded (NC) and repetition-based (RepCod) cooperative diversity, as 
described in the caption of Fig. [5] The separation vectors of both sources are 
shown in the legend of the figure. For clarity, only Monte Carlo simulations 
are shown. 



Monte Carlo simulations. An example is NC-1 shown in Fig. 
[9] and Fig.[l0l which confirms C4) in Section ITV-Dl The other 
results are not shown due to space limitations. 

Finally, we would like to emphasize that in the present paper 
many numerical results are given only as illustrative examples, 
which are aimed at substantiating analytical derivations and 
main findings. However, many other options are possible for 
network code design, which depend on the application-specific 
robustness to multipath fading requested by each source. In 
particular, it is interesting to comment on how close to the 
Singleton bound (full-diversity) the network codes considered 
in this section are. From Section IIV-CI and Section IIV-EI 
we know that for binary network codes full-diversity can be 
achieved by all the sources only for single-relay networks. 
However, full-diversity may be achieved by some sources. For 
example, the maximum achievable diversity order is equal to 
three in Fig. [2] six in Fig. [3] and four in Fig. |4] and Fig. [5] 
Similar comments apply to the other figures, but comments 
are here omitted for the sake of conciseness. The network 
code in Fig. |2] is an example where the Singleton bound can 
be achieved by at least one source. This network code can 
be considered to be near-optimum since the second source 
achieves, at least, second-order diversity. The network code 
in Fig. [3] is an example where the Singleton bound is not 
achieved by any sources, but they both have a quite high 
diversity order. From [37, Table I, row 7], it follows that other 
network codes are available to allow one source to achieve 
full-diversity order equal to six. However, the price to pay 
is that the diversity order of the other source is only equal 
to two. Thus, even though the network code in Fig. [3] is not 
strong enough to attain the Singleton bound for none of the 
sources, it provides a good diversity-multiplexing tradeoff for 
both sources. The network code in Fig. |4] offers a trade-off 
similar to the network code in Fig. [3] where no source achieves 
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Ji(s|c,c) = ! 7StD {exp(~4:-fs t Dd 2 St s)exp (Aj StD d 2 St s 2 )} 

(a) f+°° 

= / exp (-4C4 t s) exp (4(d 2 St s 2 ) {{I/is.d) exp (-C/t&b)] (31) 
Jo 

( b ) r„ ,9 /, M-l( c ) 



= [1 + 47 StB 4 t * (1 - *)] 1 « [475 fCS (1 - S )]-I" 6s *" 6s *l 



0(s|c, c,£m) = II ff?(s]c, c, £ m ,n q -i = 0) x n Qq (s| c, c,£ m , /x 9 _i = 1) 

8 = 1, Mg-1=0 9=1, M,-l=l 

e9(s|c,c,£ m ,/i a _i =l)=Er lexp^AH^^od^s^exp^A^T^fld^s 2 ) f E 9s t fl,Q (y 2 7s t fl 9 ) 



T ( s ) =E { 7sR ^ 7fig£ ,} { CX p(- 4min {7egfl,,7 J R,D}4, s ) ex P ( 4 { m ^he q R q ,lR qD }) 2 l R ] D d 2 Rq S 2 ^ 

r t ( s ) = E { 7sH , 7R?D } {exp(-4min{ 7e(zfl9 ,7 jR9jD }4" ) s)exp (4 (min { 7egflq , lRqD }) 2 7^4, fi2 ) <9 (y/^s t R q ) } 

(33) 



full-diversity but both have diversity order equal to three. One 
source could achieve full-diversity equal to four by using the 
network code in [37, Table I, row 2]. However, the diversity 
order of the other source is limited to two. Similar to Fig. 
12 the network code in Fig. [5] allows at least one source to 
obtain full-diversity. In conclusion, our numerical examples 
clearly show that using NC allows cooperative networks to 
achieve a wide range of diversity orders, which can be chosen 
in agreement with the robustness to multipath fading requested 
by each source. Furthermore, these diversity orders can be 
achieved with the same number of time-slots, and, thus, the 
same network rate. 

VI. Conclusion 

In this paper, we have contributed to the theoretical under- 
standing of network-coded cooperative diversity protocols. We 
have shown that these networks generalize repetition-based 
cooperative relaying protocols, and that they can offer, by 
properly choosing the encoding vectors at the relays, good 
design flexibility to accommodate various performance and 
rate tradeoffs. For a fixed network topology and desired end- 
to-end diversity order, the encoding vectors can be constructed 
from linear UEP block codes. By assuming the C-MRC at the 
destination, an asymptotically-tight analytical framework for 
arbitrary network topologies and binary encoding vectors is 
developed, and its achievable diversity is studied analytically. 
It is shown that the diversity order of each source coincides 
with the separation vector of the distributed network code. The 
framework accounts for two classes of relays, i.e., partial- 
and full-cooperative, and it is shown that the former class 
contributes to neither the diversity order nor to the diversity 
gain of the sources. Thus, partial-cooperative relays are of no 
use for binary network-coded cooperative diversity protocols. 
Ongoing research is concerned with the analysis of these net- 
works with non-binary modulations and non-binary encoding 



vectors, as well as with the end-to-end code design over 
more general fading channels (e.g., block-fading channels) 
with channel coding [35], [48]. 

Appendix I 

High-SNR Computation of T t (-| •) and Q ( -| •) in ([T"4l 
For high-SNR, Ft ( -| •) in (fT4l can be computed in closed- 

. , (a) 

form as shown in (13 1) at the top of this page, where: i) in = , 
it is taken into account that ~fs t D is an exponential RV with 

(b) 

parameter 7s t H ; ii) = follows from [49, Eq. (3.310)]; and 

(c) 

iii) = holds for high-SNR, i.e., 7s t D 3> 1. This concludes 
the proof of ( D31 l. 

In (TBI , let us consider the generic event £ m for 
< m < 2 Nr - 1. The cases rn = and m = 2^" - 1 
can be obtained with similar analytical steps. Then, G(-\-) 
in (TT4l can be re-written as shown in (l32l at the top of this 
page, where: i) is a short-hand to denote all channels 

from the N$ sources to R q ; and ii) it is taken into account 
that = d? R if fi q -i — and = —d 2 R if fi q -i = 1. 
To compute d32l i in closed-form, the expectations T (■) and 
T (s) = Ylt^i 9s t R r* (s) shown in (|33T > at the top of this 
page need to be studied. When computing T t (•), we assume 
9s t R q = 1 as only in this case it contributes to T (•). 

T (•) can be computed in closed-form as T (s) = Ti (s) + 
T2 (s). Ti (•) is shown in d34l at the top of the next 

(a) 

page, where: i) in =, it is taken into account that j R d 
and jeqR are exponential RVs with parameters j R d and 

(*>) ' (c) 

leqR q \ ii) — follows from [49, Eq. (3.310)]; and iii) = 

holds for high-SNR, i.e., E m /N > 1. T 2 (•) is shown 

. . (b) (c) 

in (135b at the top of the next page, where = and = 
follow from considerations similar to (134V On the other 

(ai) (02) (ai) 

hand, = and = deserve further clarifications. In = , we 
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(a) /' + 00 fQeqRq / 9 \ / o\ 

00 = J J exp (-^iJ^Drfj^aJ ex P {^CR q Dd Rq S J [(l/7fl g £>) exp (-Cfl 9 £)/7fl s £>)] [(l/%qR q ) exp (~CeqRq /%qR g )] dC,R qD dC, eq R q 



(b) 



1 + 47H zxio s(l-s) l-7K„£>/(47ij,r)7 eg H d H s(l-s)+7 e9 fl +7h o) « [4 7fl r>s (1 



(<0 



(34) 



T 2 



SeqRq 



00 '-'y T 2 ( s ;Ce<jii<j) exp (^CeqAjdf^s) [(l/7 egH J Cxp (-Ce 9 fl 9 /T^ii, ) ] ^Ce 

( a 2) /' + 00 f /' + 00 ] / \ 

^ JO \J( ( 1 /T'«9 ) EX P (-C-B, -D/7fl,D) dCitgDj exp ^-4Ce,il g d?j 9 Sj (l/legRq) exp {-CeqRq/leqR q ) dC, eqRq (35) 



(6) r 2 ,_ , . ] -1 (c) 

= l+47e,H g rf H (7e,fl,/7fl 9 D) ~ |4 7 , 



|r(NC) ,(NC)| 



eqR. q 



T 2 (s;CegHg) = / exp 4 (£ qRq /C,R qD ) d 2 Rq s 2 (l/j RqD ) exp {-C, RqD / lR q D) dC, RqD = / g{x)f{x)dx (36) 



T t (s;0) = E { 7sR9 , 7i??D } { 0X P (-4min{ 7e9j j 9i 7 jR9l5 }d^ ) s) exp (4 (min { leqRq , 7 fl,z>}) 2 l RqD d\ q s 2 ) exp (- 7St ji, /sin 2 (0))} 



(a) 



, t exp (-4 min {73,^,7'^, lR qD } d^ exp U (min {7S t fl, , 7^a, , iRqo}) ^^4,^ I ' x !> 



( = 5 rf (s; (?) + rf > (s; 0) + rf ' (s; 0) 



lS t R q 
" sin 2 (0) 



(37) 



define f (x) = (tb^i 2 ) 1 exp [-l/(xjR qD )] , #0) = 
exp {^Q qRq d Rq s 2 x^j , and T 2 (•; •) is given in (|36T > shown at 

the top of this page. From these definitions, the high-SNR 

O2) 

approximation in = follows by noticing, with arguments 
similar to [42, Fig. 1], that the behavior of g (x) around x — » 
mainly determines the high-SNR behavior of Y 2 (•). With this 

(02) , x 

in mind, = is obtained by replacing g (•) with its first-order 
Taylor expansion, i.e., g (x — > 0) « 1. 

Using the Craig representation of the Q-function [39, Eq. 
(4.2)], we have T t (s) = £ /2 T t (s; 0) dO, where T t (•; •) is 

. . (a) 

defined in d37b at the top of this page, and: i) = is obtained 
by introducing 7 



(\t) = 

eqR q 



nun T ^ i= i,2,...,iv s {ffs T fl.,7s T i?.,}, 



which is an exponential RV with parameter = 

(Ylr^xdS^Rqls^) > and ii) - follows by decom- 
posing r t (•; •) into the summation of three terms, i.e., 
r t (1) (-;-) for min {7^,75^,7^0} = ls tRq , I? 3 (•; •) 
for min {7^,7^,7^23} = 7^, and rj a) (-;-) for 
min {lS t Rq,1^R q ,lRqD\ = lR q D- Similar to the compu- 
tation of T (•), each of the three summands in = can be 
decomposed into the summation of two terms, and each of 
them can be computed in closed-form and for high-SNR by 
resorting the same analytical steps and approximations used 
to compute Ti (•) and Y 2 (■). The details of the derivation 
are here omitted due to space limitations. In particular, our 
analysis has shown that, for high-SNR, T (2) (•; •) and v[ 3) (•; •) 
decay as (E m /N )~ 2 , while (•; •) decays as (E m / Nq)^ 1 . 



Thus, rj 1 ^ (•; •) is the dominant term for high-SNR. In con- 
clusion, we have: 



T t (s; 



{3; 9) 



lS t R, 



g (l/sin 2 (0) 



(38) 

Finally, by inserting C3HE2), <E2K and CM) in <32j, Q (-| •) 
in ( fT4l can be obtained after some algebra and by taking into 

— > 1 for high- 



account that W„ R , n 
SNR. This concludes the proof. 



1 



{±leqR q ) 



Appendix II 
High-SNR ABEP of Single-Relay Networks 

Let us consider the same system model as in [22] with 
Ns sources and one full-cooperative relay (Nr = 1) with 
encoding vector gR 1 = [1, 1, . . . , 1]. By direct inspection of 
the codebook and from SectionHV] we have G^ St) = SV (St) = 
2 for t = 1,2,..., Ns- Since there is just one relay, by using 
g(-\-)in 03, the APEP (0 -> c) in O simplifies as shown 
in (1391 at the top of the next page, where: i) jVfi (s|c) = 
rii= S i_^t(s|c); ii) M 2 {s\c,£ ) =&(a| c, £ ,Mo=0); and 
iii) M 3 (s\ c,£i) = Q x <s\ c,£ 1} fi = !)■ 



From (|27>, ABEP^ V ' can be explicitly written as follows: 

N s +1 

(40) 



ABEP 



(SV) 
St 



Ns + 1 

APEP I 

T = l 



;(t>r) 



where c^*' 1 "^ denotes the codeword whose entries are c t = c T = 
1, and c p = if p 7^ t and p 7^ r for p = 1,2,..., Ns + 1. 
This simplified expression of ABEP^ V ' ) originates from the 
fact that c(*' T ) are the only codewords of the codebook with 



wu (c ( *' T) ) = 2. 
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APEP (()->•£) = (2nj) 



5+joo 
5—joo 



Mi(s\c)M2(s\c,S )s- 1 ds + (2nj) 1 



5+JOO 



M x (s|c)7W 3 (s|c,f 1 )s- 1 ds (39) 



,7T/2 


i-S+joo 


(4/tt) / 


(27TJT 1 / a 


Jo 


J (5—joo 


/.tt/2 




(4/tt) / 


(27Tj)- 1 / S 


Jo 


«/ (5—joo 



r/(9 



(a) 



tt/2 



(4/tt) / (4 + sin" 2 (6»)) 1 d<9 = (5 - y/b) / 



10 



(43) 



r/(9 



r/2 



(4/tt) / [(-4 + sin" 2 (6>)) 1 + 16 sin 4 (9) (4 - sin" 2 (0)) ' 



d9 = 1 



ABEP 



(SV) 
St 



lS t D 



1 / 1 



Ns 1 

£ 



16 \ 7i?iD T ££ 1 7.sf 



45 + V5 1 



160 



(44) 



2"R-1 / r S+joo 

APEP(O^c) = ^ (2 7 r J r 1 / 

m=0 V J*-joo 



A^i(s|c)7W 2 (s|c,f m )7W 3 (s|c,£' m )s L ds 



(45) 



By direct inspection, from ([l5T>-([T7l> we can readily obtain: This concludes the proof. 



Mi (s|£(*' r )) = {-is t Dls T D) 1 
M 2 (s\ So) =1 
Af 3 (s| £1) = (4 7e9fll )- 1 

if 1 < r < A s , and: 



[4s (1 



(41) 



Mi (s| <• 

s 



? (t,JV s +l)\ _ 



47s t DS (1 - s) 



>f 2 H Vo J =(47^*) 

+ (4 7 i?,D S (l- S ))~ 1 
1 



(4/TT) 



tt/2 



(6 1 ) + 4a) 



d9 



M : 



( S |c^+ 1 ),£ 1 )=(4 7e ^ 1 )- 1 



tt/2 



(4/tt) / (sin" 2 (<9) - 4s) d0 

(42) 
if r = A s + 1. 

By inserting fiTT i and (l42l in (|39) , the APEPs can be com- 
puted in closed-form. Due to space limitations, the details of 
the derivation are omitted. However, the adopted methodology 
is as follows: i) contour integrals are solved by using the 
residues theorem in [40, Eq. (6)]; ii) the integrals in (I4H 
and (l42l are computed only after solving the contour integrals 
in d39) , i.e., the order of integration is swapped; and iii) the 
identities in (l43l shown at the top of this page are used, where 

(a) (6) 

the equalities in = and = follow from the residues theorem 
in [40, Eq. (6)]. 

From d43) and some algebra, (|40) simplifies to (PHt shown 
at the top of this page, which coincides with [22, Eq. (32)]. 



Appendix III 

High-SNR ABEP of Partial-Cooperative Relays 



By using Q{-\-) in (02), the APEP (0 ->■ c) in (O 
can be re-written as shown in (|43T > at the top of 
this page, where: i) Mi{s\c) = nSt Ft ( s l n ) 

M 2 (s\c,£ m ) = nq=i,» q _ 1 = G<i( s \ci £ ™,fJ>q-i=0); and 
iii) M s (s\c,E m ) =Y[q=i^ q _ x =iG q {s\c,£ m ,iiq-i = !)■ 
From (05), APEP (0 ->■ cW) in (123) can be computed as: 



APEP 



(° 



(*) 



6+ j 00 



Vs 



(46) 



JS-joo ^ 1 

(b) rS+joo 

W(2 7 ri)- 1 / [47 StD s 2 (1 - s)] l ds 

(a) 

where: i) ps follows by noticing that for high-SNR the domi- 
nant addend in (|43T > is to = 0, and that M2 (s\c, £ m ) = 1 for 

every m since ' m ^ (c^) = for every m; ii) « follows 
from M^slcW) =n^ S i^(s|£ (t) ) ~ [4 7StD s(l-s)]- 1 

(c) 

with the last approximation coming from (115V . and iii) « is 
obtained by applying the residues theorem [40, Eq. (6)]. 

Let q* be the partial-cooperative relay of interest. From 
(03), APEP (0 -> c^*)) in (f23) can be computed as shown 
in ( 147) at the top of the next page because for high-SNR 
there are only two dominant addends in (143V i) m = 0; and 



IEEE TRANSACTIONS ON WIRELESS COMMUNICATIONS 



20 



APEP | 



(o^c(' 3 *)) ^(2nj)- 1 +J Mi (s|c( 9 * ),£(,) s~ 1 ds + {2TTj)- 1 +J M 2 ( s\ c( 9 *) , £ m .) Ms ( s\ c( q * ) , £ m » ) s" 1 ^ (47) 



«-t 

<5 — joo 



/ v / \ rS+joo 

(0 £<«*>) ^ (27TJ)- 1 / [^ Rq ,DS 2 (1 - »)] 

v 7 JS-ioo 



APEP 



(*>) 



(27Ti)" 1 / (47e 9 R 9 ,) _: 

*/ 5—jOQ 



ds 



7T/2 



(4/tt) / (sin" 2 (6») - 4s) l d0 



s 



(49) 



ii) m = m* with £ m * being the event that R q - is the only 
relay that forwards a wrong estimate of its network-coded bit 
to D. 



From (fTol i. we have A^2 (s| c 



of: 



.M 3 (s|c^,£„ 

x (4/tt) 



(4/tt) 



tt/2 



.) = 1 and 

+ (47fl 
r 2 (0) 



4s 



d<9 



/() 



W 2 -1 

/ (sin -2 (0) - 4s) d<2 
Jo 



Then, 
APEPfO 



(48) 

by inserting ( |48l l in d47l i we obtain 
— ► c^ 9 )) shown in d49l ) at the top of this page, 



(a) 

where: i) = takes into account that when applying the residues 
theorem only the addends containing positive poles have a 
non-zero contribution, while the others can be neglected 

0) 

[40, Eq. (6)]; and ii) = can be obtained by solving the first 
integral as in ( |46] i, and the second integral by solving first the 
contour integral. In particular, we have used the identities: 1) 
Virjy 1 //+/^ (sin" 2 (9) - 4s)" 1 s^ds = sin 2 (0), which 
originates from the application of the residues theorem [40 
Eq. (6)]; and 2) (4/tt) fi r 
(|49l , the proof is complete 



Eq. (6)]; and 2) (4/tt) f* /2 sin 2 (9) dO = 1. From © and 
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